Signatures of Thermal Dilepton Radiation at RHIC 

Ralf Rapp 

Department of Physics and Astronomy, State University of New York, Stony Brook, NY 11794-3800, U.S.A. 

(February 1, 2008) 

The properties of thermal dilepton production from heavy-ion collisions in the RHIC energy 
regime are evaluated for invariant masses ranging from 0.5 to 3 GeV. Using an expanding thermal 
fireball to model the evolution through both quark-gluon and hadronic phases various features of 
the spectra are addressed. In the low-mass region, due to an expected large background, the focus is 
on possible medium modifications of the narrow resonance structures from uo and (j> mesons, whereas 
in the intermediate-mass region the old idea of identifying QGP radiation is reiterated including 
effects of chemical under-saturation in the early stages of central Au+Au collisions. 



I. INTRODUCTION 



The identification of experimental signals from potential phase transitions in strong interaction matter associated 
with deconfinement and/or chiral symmetry restoration is the major goal of the (ultra-) relativistic heavy-ion program. 
Convincing evidence is likely to be achieved through a consistent combination of various suggested signatures. Among 
these is the copious production of thermal lepton pairs as originating from a supposedly thermalized fireball. Their 
relevance comes in two facets depending on invariant mass. 

In the low-mass region (LMR, below ~1 GeV) thermal radiation is governed by (in-medium) decays of the light 
vector mesons p, to and (p. With the dynamics in the light quark sector (u, d, s) being driven by chiral symmetry 
and its spontaneous breaking in vacuum, medium modifications of the light vector mesons should encode information 
on (the approach to) chiral restoration. The challenge is then to discriminate the thermal contribution from the 
numerous 'background' mostly consisting of meson decays after freezeout of the interacting (hadronic) system. 

At sufficiently high masses (above 3 GeV or so), primordial processes such as Drell-Yan annihilation (characterized 
by a much harder slope than typical thermal sources) or correlated decays of open anti/-charm (bottom) pairs are 
expected to prevail the spectra. Thus, the focus is on the narrow peaks from heavy vector mesons J/ V I / , T, . . ., 
in connection with deconfinement signatures through their early dissolution. Descending from high masses thermal 
radiation is believed to become competitive below the J/'F resonance. In the so-called intermediate- mass region (IMR), 
extending down to the <f>, the emission is continuum-like with a magnitude approximately given by perturbative qq 
annihilation. An enhanced dilepton yield in the IMR has been among the earliest suggestions for a signal of 
QGP formation. The idea resides on the fact that at masses substantially larger than typical temperatures, thermal 
occupation factors strongly favor contributions from the hottest (early) stages of central collisions where deconfined, 
perturbative matter ought to be formed. 

Dilepton measurements performed within the heavy-ion program at the CERN-SpS over the last decade have 
demonstrated that the above issues can be thoroughly addressed. For Sulfur and Lead projectiles both the LMR 
and IMR clearly exhibit appreciable signals above the aforementioned 'background' sources, i.e., final state hadron 
decays as well as Drell-Yan and open charm contributions, respectively. At the same time extensive theoretical analyses 
have shown that the experimental findings can be consistently understood in terms of thermal emission (coupled with 
strong medium effects in the LMR, see, e.g., ref. for a recent review), indicating that one is indeed producing QCD 
matter in the vicinity of the conjectured phase boundary. Nevertheless, the major part of the enhancement is still 
associated with the hot and dense hadronic phases in the course of the collisions. 

At the now operational collider experiments at RHIC, however, the plasma phase is anticipated to occupy much 
larger space-time volumes together with higher initial temperatures. At the same time, primordial and final hadron- 
decay sources increase. For the IMR various interplays have been elaborated in refs. and for the LMR 

assessments have been made, e.g., in ref. Jl6| . The purpose of the present article is to provide an overall picture as 
well as an improved estimate of the thermal component in both the IMR and LMR, thereby utilizing our understanding 
as inferred from the investigations performed in the context of the SpS results. 

The article is organized as follows. In Sect. [Il] we will recall basic ingredients for the description of equilibrium 
dilepton production rates from threshold to about 3 GeV. The emphasis will be on in-medium effects in the LMR, which 
will be assessed in the spectral function approach within hot hadronic matter as appropriate for RHIC conditions. 
In Sect. Ill we apply the emission rates using an expanding fireball model around midrapidity to calculate dilepton 
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invariant mass spectra, suitable for comparison with future measurements of the PHENIX detector. In the LMR 
special attention is paid to the manifestation of medium modifications of the light vector-meson resonances. In 
the IMR the predictions for thermal yields will be confronted with Drell-Yan pairs to evaluate their significance for 
identifying signals from an equilibrated partonic system. In particular, the consequences of chemical under-saturation 
in an early formed 'Gluon-Quark' Plasma will be assessed. In sect. |W we summarize and conclude. 



II. ELECTROMAGNETIC CORRELATION FUNCTION AND DILEPTON PRODUCTION RATES 



A. Definition and Spectral Decomposition 



The general object for studying thermal dilepton emission is the electromagnetic (e.m.) current-current correla- 
tor |]l7| , |l8| . In vacuum it is defined by the expectation value of the time-ordered product of e.m. currents, 

Kr(q) = -if d 4 x e**° (0|Tj eTO (z)j eTO (0)|0>. (1) 

The correlation function provides a common framework not only for different hadronic approaches but also for quark- 
gluon based calculations. This feature is particularly evident in free space where the following decomposition of its 
imaginary part can be given: 

E ((™ { v) 2 /9v) 2 lmD^(M, q; fi B ,T) , M < M dual 
. V= p,u),4> 

imnz = { (2) 

(-g^ + q^/M 2 ) (M 2 /12ir) N c £ {e q f , M > M dual . 

q=u,d,s 

From total e + e~ — > hadrons cross section data it is well established that the low-mass part of the correlator can be 
accurately described by the light vector mesons within the Vector Dominance Model (VDM), whereas above a 'duality 
threshold' of about Mdual — 1.5 GeV free quark-antiquark states satisfactorily saturate the strength (within 20% or 
so). 

The thermal dilepton production rate per unit 4- volume and 4-momentum is directly related to the (spin-averaged) 
imaginary part of the in-medium (retarded) e.m. correlator through 

= - 3 A/2 S B {qo\T) lmU cm (q Q ,q;T,fi B ) , (3) 

where the lepton tensor L^ v encodes the propagation of a virtual photon and its subsequent decay into a pair of 
leptons. 



B. Medium Effects in the Low-Mass Region (LMR) 



In the LMR medium effects arise from changes of the vector-meson spectral functions in hot hadronic matter. 
At RHIC energies, due to the low stopping of the interpenetrating nuclei, one would not expect baryons to play a 
significant role at central rapidities. However, even for a vanishing net baryon density this may not necessarily be true 
once accounting for thermal excitations of NN pairs in connection with the rather strong meson-baryon interactions 
(as compared to the relatively weaker meson-induced interactions). In the following we will employ standard many- 
body techniques to implement (hadronic) selfenergy corrections into the propagation of the three light vector mesons 
p, uj and 4>. 



1. p-Meson 



For the p-meson propagator, 



D P (M, q; T) 



M 



,(o)' 



-•pB 



(4) 
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an extensive literature about hadronic in-medium modifications has emerged over the last 10 years (see |5) for a 
recent review). Here, we incorporate finite temperature effects along the lines of ref. where the total in-medium 
selfenergy has been obtained in terms of polarizations of its pion cloud as well as resonant interactions with surrounding 
pions, kaons and rho mesons, M — n, K, K, p. In addition, modifications due to nucleons as well as the most abundant 
baryonic resonances (B = N, A(1232), iV(1440), iV(1520)) are included following refs. Since the p°-meson is 

an eigenstate under antiparticle conjugation, it equally interacts with antibaryons, thereby effectively doubling the 
impact of the baryonic component for ps close to zero. E.g., at (T, /zb) = (180,0) MeV one has qb = Qb — 0.46go 
which is appreciable (including all nonstrange resonances up to ms = 1.7 GeV as well as the lightest hyperons; 
go = 0.16 fm~ 3 denotes ground state nuclear matter density). 



2. w — M eson 



The same approach as for the p is now adopted for the us meson (see also refs. Jl9|-^T| for somewhat different 
treatments). First, its vacuum selfenergy is constructed from a combination of the direct to — > pit and the to — > 3tt 
decays derived from the interaction vertices (isospin structures suppressed) p5|,|2^] 



Cwpn = Gu>pir E/x^ctt ^ d lJ p a d T ir (5) 
Cw&n = G^ztt tp.vor ^ d v ir d a ir d T ix , (6) 

respectively. G p7rw is taken from ref. j23] which (together with a hadronic formfactor) gives the correct radiative decay 
width T u ^t T1 (using VDM) and about 50% of the full hadronic width. is then adjusted to reproduce the missing 
part. Medium effects in the 7rp cloud are introduced through thermal Bose enhancement factors (evaluated in the 
Matsubara framework) as well as the in- medium p spectral function from above, i.e., 

r w ^( s ) = i^ f A P (M) [1 + f*(w„) + f p {u p )} 2ql m i^ p7r (g cm ) 2 (7) 
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i + r{toi)} + /"K){fK) + r(to 2 )} + rmrh) 



2 2 1 / 2 2 2\ 

P1P2 - j(Po -Px -P2) 



Fcj3k(P0,P1,P2) 2 (8) 



with M m ax = \fs — win, u>i, ma x — \fs — 2m v and u>2,max — \fs — (toi + m^) (dynamical modifications of the pion 
propagators are smaller and have been neglected). The u> can also be absorbed through the inelastic reaction ton — ► -kit. 
The corresponding collisional broadening has been calculated in ref. EG] which we include via 



r w ,™ ~ (7MeV) 1-1 (9) 

(To = 150 MeV). We furthermore compute an w-selfenergy from scattering off thermal pions inducing &i(1235) 
resonance formation |^,^], which has the largest known coupling to urn states. Using the interaction vertex 

£u)7rbi = G W7rbl 6 M {g^uP ■ q ~ %pv) u v ir , (10) 

the selfenergy takes the form 

S W7r6l (M, q) = G 2 ^ [ ,f^l dX , > [fkW) - f" 1 (qo + w,(p)] F^ bl (q c m) 2 D bl (s) v^ bl (p, q) , (11) 
J (2tt) 2 2w w (p) 

where the vertex function v^b! (p, q) and the hadronic (dipole) formfactor F^bx are in close analogy to the pirai 
case [ p2[ . The coupling constant and cutoff parameter are fitted to the empirical hadronic and radiative decay widths 

In nuclear matter, according to the coupled channel analysis for vector meson- nucleon scattering of ref. fl30| ], the 
dominant coupling of the us proceeds through resonant iV(1520)iV~ 1 excitations (with a strength similar to that for the 
p meson), as well as N(165Q)N~ 1 ones. Modifications induced in the pion cloud are presumably less significant |si|, 
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especially when using soft formfactors for the irNN and irNA as suggested by the analysis of photoabsorption spectra 
on the nucleon and nuclei §53]. In line with the arguments given above (see also ref. ^4j ), we evaluate the N(1520)N~ 1 
and N(1650)N~ 1 loops at an effective nucleon density 

Qeff = Qn + Qn + 0- 5 (Qb + Qb) ( 12 ) 
to account for the effects of higher resonances as well as antibaryons. The in-medium to propagator finally reads 

Du, = [m 2 - ml + im w (r 37r + r p7r + r^^rnr) - s W7r6l - £ wB ] 1 ■ (13) 



3. 4> -Meson 



For the <fi meson various evaluations of medium effects based on hadronic rescatterings at finite temperature have 
indicated quite moderate changes of both its mass and width p7| , p8| , po| , p9[ . This also holds true for direct (pN 
interactions as indicated by photonuclear production of <f> mesons as well as the absence of baryonic resonances with 
explicitly measured <j>N decay channels^]. Therefore we do not attempt a microscopic description of this part of the <p 
selfenergy but merely parameterize the finite-temperature collisional broadening guided by the results of ref. |2Q] as 



Im S 



(T) = -m, 17" (T ) (^j (14) 

with T {) = 150 MeV, Tf ll (T ) = 9 MeV. 

In addition, there can arise width modifications from the dressing of the anti-/kaon cloud through in-medium 
<j> — > KK decays. Here even small effects in the kaon propagators can have appreciable impact on the <j) spectral 
function due to the proximity of the 4> mass to the two-kaon threshold (the p-wave nature of the <fi — > KK decay 
further enhances this sensitivity). We assess this by evaluating the 4> — > KK width in terms of in- medium modified 
(anti-) kaon spectral functions, i.e., 

M 

Im £^(M;T) = -m* T™ kR J ^4(mi) 

:(™ 2 ) fe) F^ipKf [l + /*(u;i) + /*M] (15) 



o 

M—mi 



^ v ^Zkk = ^•^ MeV: on-shell width in free space), and again neglect possible modifications of the real part. The 
(anti-) kaon spectral functions are given by 

A K M m ) = ~ 2 ^DkM™) = - 2 Im — 5^ • ( 16 ) 

m — m K - z^ K K 

For the finite-temperature part of S^-^-, which does not distinguish between kaons and antikaons, we approximate 
the results of ref. 119] as 

ImZ K<R (T) ~ -m K T K<R (T) = -m K (20MeV) 

ReE K<R (T) ~ -m K (lOMeV) . (17) 

In nuclear matter kaons and antikaons behave quite differe ntly. The KN interaction is rather weak and slightly 
repulsive. A corresponding selfenergy has been extracted as [j33 35 



1 This behavior is quite contrary to the light quark states p and ui. The reason for that is presumably related to the OZI rule: 
whereas in a p-/u>-N interaction the antiquark in the vector meson can easily annihilate and convert into excitation energy of 
a baryonic resonances, the equivalent mechanism for the antistrange quark in a <f>-N interaction is strongly suppressed. 
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T, K (g N ) ~ m K {0.13m K ) f — 



(18) 



The antikaon, on the other hand, couples to a number of hyperon excitations. A detailed recent coupled channel 
calculation |54| predicts an appreciable K width of up to 100 MeV already at nuclear saturation density, together 
with significant attraction as well. Accordingly, we employ 

T, R (m;g N ) ~ m K f^j (-60, -i 90)MeV . (19) 

Notice that the situation reverses for the interaction with antinucleons, i.e., 

Z k (qk) = X r (q n ) (20) 
Ek(qn) = ^k(Qn) (21) 

Thus, for central rapidities at RHIC, where the net baryon density in the hadronic phase is small, the combined 
contributions to the real part of both the anti-/kaon selfenergies practically vanish. This is not the case for the 
imaginary parts which will be explicitly included. 

Finally, we take into account the 13% branching ratio of the free <fi decay into irp, and treat medium modifications 
thereof in analogy to the ui — > pn decay, eq. ([?]). 

Combining all contributions yields an in- medium <j> propagator of the form 

D = [M 2 -ml- ~ £*p* - ~* ■ (22) 



4- Spectral Functions 

Fig. |] summarizes the appearance of the light vector-meson spectral functions (=imaginary part of the propagators) 
that will be employed in the calculation of dilepton spectra below: both light quark states p and u> undergo a strong 
broadening towards the expected phase boundary (one also observes a slight upward mass shift caused by repulsive 



parts in the selfenergy which are mainly due to baryonic particle-hole excitations |24 30 1 , as well as medium effects in 
the pion tadpole graph in case of the p p6| ) . The (f> seems to retain more of its resonance structure, although, at the 
highest temperature, the hadronic rescattering has increased its vacuum width by over a factor of 7 to ~ 32 MeV. 




FIG. 1. Spectral functions of the light vector mesons p (left panel) and u, (f> (right panel) in vacuum (solid lines) as well 
as in hot baryon-poor hadronic matter as expected under RHIC conditions: (T, /ijv) = (120,91) MeV (long-dashed lines), 
(T,un) = (150,40) MeV (dashed-dotted lines) and (T,fj, N ) = (180,27) MeV (short-dashed lines). For the definition of g e}} see 
eq. @. 
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It is interesting to note that at T = 180 MeV, where for the chosen value of ps = 27 MeV the net baryon density 
amounts to only 0.14po, the combined effect of baryons and antibaryons causes an appreciable further broadening of 
both the p and uj spectral functions as compared to the pure meson gas results, see fig. ||. This does not apply for 
the 4> meson. It would be very illuminating to scrutinize these features with high resolution/statistics experiments in 
more baryon-dominated systems as created in fixed target heavy-ion collisions. 




FIG. 2. p (left panel) and ui, (j> (right panel) spectral functions at T — 180 MeV with (short-dashed lines, ps = 27 MeV) 
and without (dashed-double-dotted lines) the effects of anti-/baryons as expected for central rapidities at RHIC. 



C. Intermediate Mass Region (IMR) 

The nonpcrturbative mechanisms related to the deconfinement/chiral restoration transition necessarily imply sub- 
stantial medium effects for the light hadronic bound states, as discussed in the previous section. Beyond invariant 
masses of M ~ 1.5 GeV, however, the interactions in the free vector correlator essentially become perturbative, cf. 
eq. (||). Naively, one might anticipate in-medium corrections to be of order a s or 0(T/M), 0(p q /M), but theo- 
retically the impact of 'soft' corrections to the finite-temperature qq rates is not very well under control, see, e.g., 
refs. [39| - |42[ . For the (closely related) Drell-Yan annihilation process, corrections of similar nature are subsumed in 
the famous ii"-factors, which, for modern parton distribution functions (e.g., GRV-94 p3|), are around 1.3-1.5. This 
should roughly reflect the uncertainty when employing the lowest-order thermal qq — > Fl~ annihilation for the QGP 
phase as done below (and might be considered as a lower limit of the true emissivity). Making furthermore use of 
'duality' arguments p7|J^j5^] we will in fact apply the corresponding dilepton production rate Jig] , 

d s N^ _ a 2 T B ^ 2 (s_ + y) {x+ + exp[-p q /T]) 

also in the hadronic phase, extrapolating it down to masses of about 1 . 1 GeV. The latter is motivated by the observation 
that lowest-order in temperature effects |46[ already amount to a reduction of the 'duality threshold' towards the 
resonance (see, e.g., the discussion in 




D. Thermal Rates 

In fig. U we compare the 3-momentum integrated dilepton production rates in thermal equilibrium, 

dR ee _ f d 3 q dR ee 

dM 2 ~ J 2q d 4 q 1 ' 

at fixed temperature and density for 3 different ways of evaluating the electromagnetic correlator: using hadronic 
degrees of freedom with (solid line) or without (dashed lines) medium effects in the vector-meson spectral functions, 
as well as perturbative quark- antiquark annihilation (dashed-dotted lines). 



G 




FIG. 3. 3-momentum integrated thermal dilepton production rates as appropriate for the hadronic phase at RHIC (left 
panel: (T, /ijv) = (150,60) MeV, right panel: (T, fir?) = (180,27) MeV). The dashed lines are obtained when using the free p, 
ui and (j) spectral functions, whereas the full lines follow upon inclusion of hadronic medium effects as described in the previous 
section; the dashed-dotted is the result for lowest order q s quark-antiquark annihilation. 

The hadronic calculations have built in contributions from 2- and 3-pion as well as K K states through the p and 
u> as well as the cj>, respectively. These are sufficient to saturate the strength in the electromagnetic correlator from 
the two-pion threshold up to about 1 GeV. Beyond, 4-pion states start to become dominant, most notably through 
■na\ and iru) annihilation pa), which are rather structureless and merge into the aforementioned quark-hadron duality 
from Mduai > 1.5 GeV on(in vacuum). 

The medium effects in the spectral functions are directly reflected in the dilepton rates. The thermal Bose factor 
in the latter augments the enhancement on the low-mass side of the resonances (a factor of 3-10 over the free case 
around invariant masses M ~ 0.4 GeV for T = 150 — 180 MeV). To a large extent, the underlying processes can be 
associated with (in matter) Dalitz decays u), ai(1260) — > ne + e~ or N* — > Ne + e~ , being proportional to the imaginary 
parts of the vector-meson selfenergies. At the same time, a strong broadening of the resonance structures is implied 
as a result of resumming the large imaginary parts [^J2^| in the propagators. Taken together, both features seem to 
drive the hadronic in-medium curve towards perturbative qq annihilation: at T — 180 MeV the magnitude and slope 
of the former in the M ~ 0.5 GeV region are quite reminiscent to the plasma rate (not so at T = 150 MeV); even in 
the to resonance region, where the maximum of the free rate is suppressed by a factor of ~6, the in-medium hadronic 
and the qq result agree within a factor of 2. On the other hand, at T = 150 MeV where the in-medium suppression 
amounts to a factor of ~3, in-medium hadronic and qq rates differ by a factor of ~4. 

One concludes that the indications for quark-hadron duality in dilepton production close to the phase boundary, 
as put forward in ref. Q for SpS conditions, also (approximately) persist in the (almost) net baryon-free regime at 
RHIC. 



E. 'Background' Sources 



To be able to extract information on thermal dilepton radiation in the complicated environment of heavy-ion 
reactions a reliable assessment of the multiple other sources is mandatory. At low masses these are mostly Dalitz 
decays of 7r°, r\ and u> mesons as well as the direct decays of produced uj's and tfi's in the final state (commonly denoted 
as 'hadronic cocktail'). Restricting ourselves to invariant masses above 0.5 GeV, we explicitly compute the latter two 
contributions based on the lo- and ^-abundances as obtained from the thermal fireball model at thermal break-up, 

^gp = Fv{M)2 V '° J ^ fV ^ T ^ ' (25) 

Here, Vf is the (final) fireball 3-volume, f v are Bose distributions and 
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F v {Mf = m% \D° V {M)\ 2 (26) 

the (free) electromagnetic formfactors (V = u>,<f>). Note that in such a framework the contribution from p-mesons, 
due to their short lifetime (r p < t™ c = 1.3 frn/c), is more appropriately evaluated through a slight extension of the 
fireball evolution, since their medium modifications around frcczeout will still be significant. The same argument 
holds for higher resonances in the M ~ 1.5 GeV region (ai, p' , etc.), which effectively make up the dual rate, eq. (|2^). 

For still higher masses (above M > 1.5 GeV) the most important hadronic contribution arises from the associated 
semileptonic decay of charmed mesons (D, D — > e + vX 1 e~vX) which originate from primordial production of cc pairs. 
A comprehensive evaluation of its features has been undertaken in ref. || , as well as in the follow-up studies in pl| , ^2[ . 
We will therefore not attempt to recalculate it but include the results and arguments given in these references in the 
discussion of the spectra below. 

Finally, at masses starting from 1 GeV (primordial) Drell-Yan (DY) annihilation constitutes an important source 
of dileptons. Again, a detailed computation of DY-pairs at RHIC energies has been performed in ref. ||. We here 
restrict ourselves to the case of central collisions of two equal nuclei with mass number A, 

riN AA 3 rlrr NN 

DY (b = 0)= A 4/3 DY (27) 

dMdy [ ~ ' ~ 47ri?2 A dMdy ' [ 

in terms of the (standard) elementary Drell-Yan cross section in a nucleon-nucleon collision, 

5jS = a 'Sf £ 4 hew*) + few*)] • <*> 

q— u,d,s 

The (anti-) quark distribution functions (7(21,2) (9(21,2)) are taken from the leading order (LO) GRV-94 parameteri- 
zation |43f| (including explicit isospin asymmetries) in connection with a if- factor of 1.5 extracted from p-A data [Ejj . 
Their arguments are related to the center-of-mass rapidity y and invariant mass of the lepton pair as 21,2 = 2e ±y with 
2 = M/y/s, where s denotes the total cms energy of the nucleon-nucleon collision. In e q. (p7| ), the root-mean-square 
radius parameter i?o — 1-05 fm arises from a folding over a Gaussian thickness function |48[. For (slightly) off-central 
collisions as considered below we simply replace A in eq. ( p7| ) by the number of participants. Effects of nuclear 
shadowing have been estimated to reduce the total DY-yield in the RHIC regime by 20-50% ||, which is within a 
similar uncertainty as implicit in the thermal production rate, eq. (p3h. 



III. DILEPTON PRODUCTION FOR Af=0.5-3 GEV AT RHIC 



A. Thermal Fireball Model 



The total dilepton yield from thermal radiation in a heavy-ion collision obviously depends on the space-time history 
of the system, which the thermal rates have to be integrated over. Here we employ a simplified thermal fireball 
description including, however, the most important features of a full hydrodynamic simulation in a transparent way. 
Assuming a spatially averaged (isotropic) temperature-density profile at each time incident, the time- and momentum- 
integrated thermal mass spectrum can be written as 

tf a 

dN th r P Afd 3 n dR 

^f = j dt V Fc{t) J ^-[M,q;T(t),^(t),ti B (t)} Acc(M, q; y ) . (29) 



Experimental acceptance corrections are encoded in the factor Acc(M, g;y ), which is prescribed by the specific 
detector characteristics, and depends on the kinematics of the produced pair as well as the rapidity positioning 
of the thermal fireball (yo = at RHIC). The latter has to be chosen in accord with measured hadronic rapidity 
distributions. As thermal fireballs typically imply smaller spreads in particle rapidity than found in heavy-ion reactions 
at ultrarelativistic energies Q), multiple fireballs have to be used. At the full SpS energy, e.g., the full-width- 
half-maximum of negatives, dN h - /dy, amounts to about r^ pS ~ 4. For thermal pion distributions at the relevant 
temperatures one has T^ 1 ~ 1.5 — 2, so that two fireballs are minimally required for SpS conditions, see, e.g., ref. p4[ . 
At RHIC energies, a larger number would be necessary for a complete description. However, since our focus here is 
on a comparably small window around midrapidity (Ay = ±0.35 as covered by the two electron arms of the PHENIX 
detector), a single fireball suffices. The time-dependent volume is modeled by a cylindrical expansion as 
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V FC (t) = (z + v z t) tt (r + 2 a ^ <2 ) 2 



(30) 



where r Q is the initial transverse overlap of the two colliding nuclei (determined by the impact parameter b), v z — c 
the longitudinal expansion velocity (reflecting the approximate rapidity coverage), and a± = 0.03 fm 2 /c 3 is chosen to 
give a final transverse flow velocity of v± ~ 0.6c in connection with a typical freezeout time of tf a ~ 20 fm/c. The 
largest sensitivity of our results is attached to the initial longitudinal size zq being equivalent to a formation time (or 
initial temperature), which will be explicitly addressed below. 



1. Hadronic Phase 



To determine the thermodynamics of the fireball we employ the picture of subsequent chemical and thermal freeze- 
out, which has proven very successful in describing many observables in heavy-ion collisions over a wide range of 
energies (from SIS to SpS), such as particle multiplicities |50|]5l| l, flow, HBT radii, etc.. From the chemical freezeout 
point on the composition of hadrons which are stable on the time scale of the fireball lifetime does no longer change. 
In systems with relatively low baryon content (SpS and RHIC) it practically coincides with the critical temperature 
of the chiral/deconfmement transition. Thus, for RHIC we fix it at T c h = T c ~ 180 MeV. The two necessary in- 
put numbers are then the rapidity densities of charged particles and net baryons; dynamical models suggest, e.g., 
dN ch /dy ~ 800 and dN% et /dy ~ 20 at midrapidity for the full RHIC energy of = 200 AGeV @ (the charged 
particle multiplicity is in line with a naive extrapolation based on first RHIC measurements at lower energies J52| ). 
With a standard hadronic resonance gas equation of state (EoS) these numbers imply an entropy per baryon of 
s/g 7 g et = 250, corresponding to a baryon chemical potential of fig 1 = 27 MeV at T c h- Isentropic expansion (i.e., at 
constant s/g 7 ^ 1 ) then governs the evolution of (T, /is) towards thermal (kinetic) freezeout, cf. left panel of fig. |J. 
Since pion and kaon numbers are, by definition, also conserved, pertinent (effective) chemical potentials \i v and [Lk 
are simultaneously incorporated. As opposed to SpS conditions, their values at thermal freezeout for RHIC are quite 
small reaching ^ ~ 10 MeV, ^ ~ 35 MeV at T th ~ 130 MeVQ. Nevertheless, we include their enhancement effect 
in the thermal dilepton production rate from the hadronic phase as well as in the cocktail contribution, cq. (B5H, via 
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FIG. 4. Thermal fireball model description of central heavy-ion collision at SPS and RHIC energies; left panel: Thermo- 
dynamic trajectories for isentropic expansion in the temperature and nucleon-chemical-potential plane with additional tt- and 
Tf-number conservation in the hadronic phases; the short-dashed lines enclose a conjectured region for the phase boundary 
between hadronic and quark-gluon matter, the long-dashed line represents an estimate for the thermal (kinetic) decoupling 
conditions in heavy-ion reactions, and the dots are the chemical freezeout points as extracted in a thermal model analysis for 
fixed target experiments at the SpS (160 AGeV) §o), AGS (11 AGeV) §(J and SIS (2 AGeV) @; right panel: time-projected 
temperature evolution from the formation time to thermal freezeout for central Au+Au (Pb+Pb) collisions at RHIC (SpS). 



The reason for that is the larger number of baryonic resonances at SpS, which serve as a storage of pion-number in the earlier 
stages, being released into pionic degrees of freedom towards thermal freezeout. 



9 



additional fugacity factors, which are z\ for p decays, z% for oj decays, z% for <f> decays, and z\ in the IMR corresponding 
to nai and ttlu annihilation [z x = exp[[i x /T}). 

We finally quote the total and net baryon densities at chemical freezeout for the fireball evolutions shown in Fig. 
at RHIC one has Qg et c h = 0-14{?o and Qg*^ — 0.54g> , whereas at SpS £>a* c/l = l-5f?o and g^lh = l-41p - 



2. Quark-Gluon Plasma and Mixed Phase 

For the early stages (starting at the formation time To) we assume an ideal QGP EoS characterized by an entropy 
density 



S QGP 



An 2 



l QG 



90 

(with a quark-gluon degeneracy cLqg = 16 + lQ.bNf and Nf — 2.5), ignoring effects from a net baryon number 
(/iq -C T). Since most of the entropy is presumably produced in the pre-equilibrium stages of a collision with little 
changes once thermal equilibrium is established (541, we can extrapolate from the chemical freezeout backward into 
the QGP phase using isentropes as before. At T c , the fraction / of matter in the hadronic phase can be inferred from 
the standard entropy balance 



(31) 



/ shg(T c ) + (1 - /) s QGP (T c ) = S tot /V FC (t) , 



(32) 



while for s(t) = Stot/VFc(t) > sqcp(T c ) the entire system is in the plasma phase. 

At collider energies (RHIC and LHC) gluons are expected to play the dominant role in early entropy production 
and equilibration processes Q , due to both their larger color charge as well as the steep increase of their structure 
functions towards low (Bjorken-) x. Although thermal equilibration seems to be achieved at a formation time of well 
below 1 fm/c, the chemical abundances of gluons and especially (anti-) quarks are likely to be strongly under-saturated. 
E.g., in the 'self-screened parton cascade' model (SSPC) the initial parton fugacities at the thermal formation time 
r = 0.25 fm/c amount to A s (t ) = n g {T )/n e g q (T ) = 0.34 and A ? ,g(T ) = n q ^{T )/n q q q {T ) = 0.064 Q (estimates 
based on perturbative QCD are even lower, see, e.g., ^7|). To study the impact of such a scenario on the dilepton 
radiation from the plasma, we employ the results of ref. KJ), where the SSPC initial conditions have been evolved 
using hydrodynamic rate equations until the onset of hadronization. For our purposes we parameterize the time 
dependence of the fugacities as found there by 

(33) 
(34) 

(t in [fm/c]), cf. fig. § 
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FIG. 5. Time evolution of anti- /quark (squares) and gluon (dots) fugacities as computed in a hydrodynamic simulation 
based on initial conditions obtained from a parton cascade model |Hq ]. The dashed lines show the results of the analytic fit, 
eqs. (Eh and (63). 



10 



Entropy production after To was found to be on the 10% level which will be neglected here. Thus we can calculate 
the temperature (needed for the dilepton production rate) from the entropy density via 



1 1 47T 2 

4gpW = I 16A s(*) + ^.5N f X q (t)} — T 3 



90 

~90" 



d QG(t) — T 3 . (35) 



Since chemical equilibration is not attained before entering the mixed phase, the construction of the latter has to be 
modified. We assume its onset to be determined by the same critical entropy density as in the equilibrium scenario 
so that the critical temperature 



T QGP {t) = 



on c c l 1 I 3 

yu b QGP 



4 d QG (t) 



(36) 



acquires a time dependence while for the hadron gas T^ G — 180 MeV is kept fixed. Eq. (|32| ) still serves to obtain the 
volume fraction of the hadron gas, fit). 

As a consequence of the quark-undersaturation, the equilibrium dilepton production rate from qq annihilation, 
eq. (p3|), is (approximately) proportional to a fugacity factor X q Xq, which translates into a strong suppression of 
the thermal yield. However, if the same initial entropy density is created as in the equilibrium scenario (which 
largely governs the finally observed number of secondaries), this suppression will be (partially) compensates by higher 
temperatures in the thermal Bose factor of the virtual photon 0. Also, if the gluon fugacities are substantially larger 
than the anti-/quark ones, a s -corrections to the qq rate become relatively more important, e.g., through the process 
qg — > qee. As shown in ref. |l5| related contributions are comparable to qq production in a more extreme scenario of 
undersaturation with X g /X q — 10. Since in our case X g /X q < 3 after r = 1 fm/c (and quickly decreasing further), we 
keep ignoring a s -corrections (as done for the equilibrium case) being aware of possibly underestimating the thermal 
yield from the off-equilibrium plasma by a few tens of percent. 

B. Dilepton Spectra 

Let us now turn to the results for the space-time integrated dilepton spectra in central Au+Au collision at y/s — 
200 AGeV. Our standard parameter set for the volume expansion within the thermal fireball model, eq. ([50]), is chosen 
as v z = c, ?'o = 6.5 fm (reflecting an impact parameter b ~ 1 fm), a± — 0.03 fm 2 /c 3 as well as initial longitudinal size 
zq = 0.6 fm (translating into an initial temperature of To = 370 MeV for a chemically equilibrated QGP) and total 
fireball lifetime of t/ = 20 fm/c (resulting in a thermal freezeout temperature of T t h — 133 MeVQ). 

We first focus on the LMR, cf. fig. ^. Only lepton pairs which fall into the rapidity range —0.35 < y < 0.35 are 
included, integrated over all transverse momenta. The main significance of in-medium effects in the radiation from the 
thermal fireball (left panel of fig. ||) are a (strong) suppression of the yield in the p-omega region and an enhancement 
below M=0.65 GeV as well as above 0.85 GeV. The 4>, on the other hand, survives as a pronounced, albeit broadened, 
resonance structure. These features become less distinct once the 'cocktail' contribution (from lu and <p decays) and 
the smooth yield from the QGP (which constitutes around ^20%) are added. Between M = 0.7 GeV and 0.8 GeV 
one still observes an appreciable ^50% reduction caused by in-mcdium effects (and a factor > 2 enhancement towards 
M = 0.5 GeV). Experimentally it might be even possible to subtract the 8 MeV wide contribution of free u> decays 
and then try to extract a remaining structure of ^50 MeV width - originating from u> decays within the interacting 
fireball - which is still well below the free p width. In the tj> region our results suggest that such a procedure is less 
promising to succeed. The factor of 2 increase just below 1 GeV might be more easy to identify. 



B note that this value is about 10-15% larger than at SpS energies owing to the substantially smaller baryon content at RHIC 
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FIG. 6. Low-mass dilepton spectra around midrapidity in b = 1 frn Au+Au collisions at RHIC energies. Left panel: thermal 
radiation from vector-meson decays in the hadronic phase with and without medium effects in the pertinent spectral functions, 
as well as from perturbative qq annihilation in the plasma phase (dashed-double-dotted line). Right panel: combined yields 
from the hadronic and plasma phase as well as from ui, (f> — * ee decays after freezeout; full and dashed-dotted lines correspond 
to the calculation with and without medium modifications, respectively, in the p, u> and 4> spectral functions. 

Next, we move to the IMR, first assuming QGP formation and evolution in full chemical equilibrium, cf. fig. [?]. 
From the left panel one sees that the hadron gas radiation prevails until M ~ 1.5 GeV before the thermal radiation 
from the QGP phase takes over due to its smaller slope parameter (originating from larger temperatures). The QGP 
signal continues to outshine the (primordial) Drell-Yan yield up to about M ~ 3 GeV (similar features have been 
found in refs. [^9|,^5|). Taken together, the total sum exhibits a noticeable change in slope between 1 and 4 GeV. 




FIG. 7. Intermediate-mass dilepton spectra at RHIC energies around midrapidity assuming a chemically equilibrated QGP 
throughout its lifetime. Left panel: Decomposition of the thermal fireball radiation (with To = 370 MeV, t/ = 20 fm/c) into 
QGP (short-dashed line) and hadron gas (long-dashed lines) parts (including their respective yields from the mixed phase), 
compared to Drell-Yan annihilation (dashed-dotted line) and the total sum (note that open charm decays are not accounted 
for). Right panel: Influence of the initial temperature and fireball lifetime on the plasma and hadron gas radiation, respectively. 
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Not included in fig. [7] are the contributions from semileptonic decays of correlated cc pairs. Extrapolations of 
their production cross sections in pp collisions to Au+Au collisions at RHIC lead to rather large dilepton emission, 
exceeding thermal radiation estimates by up to an order of mag nitude HQ at invariant masses M ~ 2 GeV. However, 
subsequent analysesj O, l^,[l4| investigated energy- loss effects for the charm quarks when propagating through hot 
and dense matter |6C|3l]. It was shown that the (rather hard) primordial transverse momentum distributions are 
substantially softened, becoming essentially thermalized. Consequently, the pertinent lepton pair invariant mass 
spectrum is strongly depleted above M ~ 2 GeV, by about two orders of magnitude, thus dropping below the Drell- 
Yan yieloQ. In addition, one might be able to assess the correlated charm decays independently through efi coincidence 
measurements |6^] or secondary vertex displacements |6^] . 

To test the sensitivity of our results with respect to the fireball evolution we have displayed in the right panel of 
fig. [?] the outcome of the following two calculations: first, doubling the initial longitudinal extent (or formation time) 
to zq = 1.2 fm/c (which translates into a reduction of the initial temperature to Tj = 295 MeV), the QGP radiation 
is appreciably decreased (by ~50%). In turn, the hadron gas evolution and corresponding dilepton yield is practically 
identical to the zq = 0.6 fm case. Second, doubling the transverse acceleration to a± = 0.06 together with a reduced 
lifetime of tf — 15 fm/c (to maintain a similar freezeout temperature of Tth — 135 MeV as before) reduces the hadron 
gas part accordingly (by ~25%), whereas the plasma yield is essentially unchanged. The latter can be attributed to 
the predominant longitudinal expansion in the early stages leaving little impact of an increased transverse acceleration 
on the first 5 fm/c or so. Thus, not surprisingly, the magnitude of the plasma yield is closely linked to the formation 
time, whereas the hadron gas yield is a rather direct measure of the lifetime the system spends in the hadronic phase. 

Furthermore, we address the question to what extent chemical undersaturation of thermal quark and gluon distri- 
bution functions in the early plasma stages affects the IMR dilepton spectra. The time evolution of the fugacities 
(cf. fig. ||) is implemented as described in the previous section (to account for the twice larger initial pressure we also 
used an accordingly increased transverse acceleration). Requiring equal initial entropy densities, the off-equilibrium 
scenario leads to a slightly reduced (enlarged) QGP yield below (above) M ~ 3 GeV due to a harder slope caused by 
the higher initial temperature^, see fig. ||. 




FIG. 8. Consequences of a chemically undersaturated QGP for IMR dilepton spectra at RHIC energies. Left panel: compar- 
ison of the plasma yield from a chemical equilibrium (short-dashed line) and off-equilibrium (full lines) scenario adopted from 
refs. (the hadron gas contribution is virtually identical for both cases). Right panel: total sum of hadron gas, Drell-Yan 

and QGP contributions with (full line) and without (dashed line) chemical undersaturation in the latter. 



4 Most of the open charm contribution is redistributed to M < 1 GeV; this reflects the kinematics of the two leptons emerging 
in randomly oriented semileptonic 3- and 4-body decays of c(c) — > I v{y)X, where the (thermal) charm quarks carry little 
kinetic energy so that Ei ~ m c /3 ~ 0.5 GeV. 

5 The small deviations from the results presented in ref. p3] arise from the fact that in there the entropy of the QGP was 
calculated with Nf — 3 (rather than Nf = 2.5 as here). 
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This result is a little less optimistic than the one of Shuryak and Xiong M], who predict an enhancement of thermal 
IMR dilepton radiation in a similar 'hot-glue' scenario (see also ref. §39]). On the other hand, Lin and Ko, using 
initial conditions estimated from perturbative-QCD |64| , find very strong suppression of the thermal yield which 
comes out well below Drell-Yan annihilation fl5|] . The employed initial conditions, X q o = A^o = 0.0072, X g o = 0.06 
and Tq = 570 MeV, correspond to an entropy density of sq ~ 12 fm~ 3 , which is an order of magnitude smaller than 
in the SSPC approach adopted here (eqs.(^), (0), and ( |35|) give s°/^ = 122 fm~ 3 ; note that for a fully equilibrated 
QGP at T = T c ~ 180 MeV one has s = 14 fm" 3 ^ 

Finally, in fig. [| we summarize the results for the combined LMR and IMR spectra with additional acceptance 
restrictions on the single-electron tracks as applicable for the PHENIX detector (pf > 0.2 GeV, \y e \ < 0.35, opening 
angle O ee > 35 mrad and an overall pair mass resolution of <tm/M = 0.5 %). 




FIG. 9. Total dilepton spectrum from b = 1 fm Au+Au collisions at RHIC energies around midrapidity including schematic 
experimental acceptance cuts appropriate for the PHENIX experiment (the w,<j> — > e + e~ cocktail contributions are not sepa- 
rately shown but included in the solid curve). Note that semileptonic decays of correlated anti-/charm and anti-/bottom quarks 
have been left out. 



IV. SUMMARY AND CONCLUSIONS 



Within a thermal framework we have discussed properties of dilepton radiation from ultrarelativistic heavy-ion 
collisions under conditions relevant for the PHENIX experiment at RHIC. 

In the first part of this article dilepton emission from static equilibrium matter has been evaluated. The emphasis was 
on the low-mass region (M < 1 GeV) , where it was found that medium modifications of the light vector mesons induce 
an appreciable reshaping of the free electromagnetic current correlator. Based on standard many-body techniques, 
hadronic interactions in hot matter at small net-baryon content generate a strong broadening of both the p and 
ui spectral functions. Close to T c , the lo width increases by almost a factor of 20 over its vacuum value, with 
non-negligible contributions from rescattering off thermally excited anti-/baryons. The <f> meson retains more of its 
resonance structure, although its free width is also augmented by up to a factor of 7-8. With that exception in 
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mind, the resulting in-medium hadronic dilepton production rates approach the continuum shape of lowest-order qq 
annihilation to within a factor of 2 in the vicinity of the phase boundary. We have thus argued that towards T c 
'quark-hadron duality' in the vector channel, which in vacuum is marked by the applicability of perturbative QCD 
above M ~ 1.5 GeV, extents down to rather low masses of ~0.5 GeV (this scale is in fact comparable to typical 
thermal energies of free (massless) partons). These findings comply with earlier analyses that were carried out for 
more baryon-rich matter. 

To calculate dilepton spectra for central Au+Au collisions at RHIC energies we have employed a thermal fireball 
model including chemical off-equilibrium effects in both the QGP and hadron gas phases. In the hadronic evolution 
from chemical to thermal freezeout quite moderate pion and kaon chemical potentials of 10-30 MeV emerge. The 
resulting low-mass dilepton yield reflects upon the lifetime of the hadronic system. As a promising signature, the 
final invariant-mass spectra contain a ~50 MeV-wide structure from in-medium us decays. This is well below the free 
p width (150 MeV), but much larger than the width associated with free us decays (8.4 MeV). The latter, leading to 
copious dilepton production after freezeout, will mask the in-medium contribution in the measured spectra. However, 
with sufficiently high mass resolution as anticipated in the PHENIX experiment, the u> 'cocktail' yield should be 
subtractable, thus enabling access to the in-medium part. 

In the intermediate-mass region the focus was on radiation from a thermalized partonic medium. Here, our in- 
vestigation was of a more qualitative nature. On the one hand, this is due to uncertainties in the QGP production 
rate (which we might have underestimated by up to a factor of 2) as well as in the formation/composition of the 
'Gluon-Quark' plasma. Nevertheless, our main point in this context should not be affected, namely that chemical 
undersaturation scenarios for gluons and especially anti- / quarks in the early stages do not decrease the thermal radi- 
ation in the 1.5-3 GeV region below the Drell-Yan yield - provided the production of the finally observed entropy (as 
encoded in produced secondaries) essentially occurs within the first fraction of a fm/c after nuclear impact (parton 
cascade models and hydrodynamic simulations support this assumption). On the other hand, the question remains 
whether the dilepton yield from correlated open charm decays can be determined accurately enough to disentangle 
QGP signals. A combination of experimental {e.g., ep coincidence measurements) and theoretical efforts {e.g., reliable 
estimates of energy loss effects in charm quark propagation) might well achieve the required sensitivity. 

Finally, let us comment on further issues that are likely to be addressed with potential upgrades of the PHENIX 
detector (such as additional vertex tracking devices [p3| ). First, a direct measurement of open charm production would 
improve the possibilities of studying direct QGP radiation substantially. Second, an evaluation of combinatorial 
background could facilitate a measurement of the low-mass continuum around the vector-meson resonances, thus 
scrutinizing in-medium modifications of the p meson in a net baryon-free environment. Since the isovector (p) channel 
has a well-defined chiral partner in form of the a\ channel, questions concerning chiral symmetry restoration can be 
investigated. Similar connections are less direct for u> or <j> mesons. 
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